Dynamic scattering channels of a double barrier structure 



M. Moskalets 1 - 2 and M. Biittiker 1 

1 Departement de Physique Theorique, Universite de Geneve, 
CH-1211 Geneve 4, Switzerland 
2 Department of Metal and Semiconductor Physics, 
National Technical University "Kharkiv Polytechnic Institute", 
61002 Kharkiv, Ukraine 

(Dated: April 2, 2008) 

We calculate analytically the Floquet scattering matrix for a periodically driven double-barrier 
structure. Our approach takes into account dynamical effects which become necessarily important 
when electrons propagate through a system subject to a fast drive. It is convenient to represent 
the Floquet scattering amplitude as a sum of amplitudes corresponding to different times spent by 
an electron inside the structure. These amplitudes define dynamical scattering channels. Then we 
represent the dc current generated by the structure as a sum of two generic contributions. The first 
one is due to photon-assisted interference processes within the same channel and the second one 
is due to inter-channel interference processes. At zero temperature both contributions are present 
while at high temperatures and/or high driving frequencies only the former survives. 

PACS numbers: 72.10.-d, 73.23.-b 



I. INTRODUCTION 

Conducting structures subject to periodically varying volt- 
ages or magnetic fields are of a fundamental interest. A 
phase-coherent conductor subject to a slow ac drive^ un- 
der quite general conditions 2 - becomes a quantum pump^ - — 
which is able to generate a dc current. Such a current gen- 
erated by an unbiased mesoscopic conductor was detected 
experimentally.— - — Theoretically the quantum pumps are 
analyzed for both slow^ 2 - - — and high2£ - — driving frequen- 
cies. 

Compared to adiabatic pumps, high-frequency pumps have 
the advantage that the currents which are generated are 
much larger. In addition, at low frequencies a proper pump 
current is difficult to distinguish from a current generated 
by rectificatio n 19 ' 21 . In contrast high-frequency pumps have 
been realized in recent experiments^ - —. These pumps might 
find an application for metrological purposes^ and or as an 
on demand electron source for solid state systems for future 
quantum information processing. Possibly another applica- 
tion is the generation of correlated particles^ - — of a pump 
outside the quantized pumping regime. 

Though several analytical results are available^ - -^, the 
theoretical analysis of pumps at high driving frequencies is 
based mainly on numerical calculations^ - —. The aim of the 
present paper is to develop an analytical approach allowing 
to describe the scattering by mesoscopic systems at arbitrary 
driving frequencies. We consider a model consisting of two 
barriers with varying strength and with a uniform varying 
potential in between them. Such a structure can serve as a 
simple model of an actual structure used in experiments^ - — 
and can serve as a prototype of a quantum pump with im- 
movable in space potentials. Our approach is based on the 
Floquet scattering matrix formalism which we developed in 
Ref.Q 

The paper is organized as follows. In Sec(TI]we represent 
the Floquet scattering matrix of a fast driven double-barrier 
system as a sum of the contributions arising from different 
reflection/transmission processes accompanying an electron 
propagation through the system, Fig.[T] Our approach takes 



(a) 

x=0 x=L 



x=0 x=L 



S ir , 11- 



1 



+ 



t 

t-2x 



1-T 



(b) 



Sin,21- ' 



x= 


x= 


L 


t-X 




t 


—> 




— a- 



+ 




+ 



t-3x 



FIG. 1: Some processes contributing to reflection amplitude 
(a) and to transmission amplitude Sin, 21 (b) are de- 
picted. The point-like scatterers located at x — and x = d 
are shown by vertical lines. The particle leaves the system 
at time moment t. We also indicate the time moments when 
the reflection/transmission at point-like scatterers occurs, r 
is the propagation time from one point-like scatterer to the 
other. 



into account effects due to finiteness of the time spent by 
electrons inside the system. In Sec lIIII we find the dc gen- 
erated current at arbitrary driving frequency and arbitrary 
temperature. The correlation properties of currents which are 
generated by the pump are explored in Sec lIVI The currents 
flowing through the structure which in addition is subject to 
an external dc or ac bias are calculated in Sec(V] The discus- 
sion of our results and a conclusion are given in Sec lVIl 



II. FLOQUET SCATTERING MATRIX 

It is convenient to introduce the matrix Si n (t,E) which 
jives the current amplitude of particles exciting the scatterer 
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at time t if they have been incident (index in) on the scatterer 
with energy E. The Fourier coefficients of this matrix define 
the Floquet scattering matrix Sf- 

T 

S F (E n ,E)= J ^e inM S in (t,E). (f) 
o 

Here T = 2ir/Q is the period of a parameter which describes 
the modulation of the scattering matrix, and E n = E + nhfl, 
is the energy of carriers which have absorbed n modulation 
quanta, n is an integer. 

For hfl <C E one can obtain an analytical expression for 
Si n (t,E) of a periodically driven system consisting of two 
point-like scatterers placed at a distance d of each other with 
a potential U in between them. For a structure with contacts 
labeled a,/3,.. the matrix element Sin,a@{t, E) can be given 
as a sum over different paths C^l which an electron incident 
in contact j3 can follow to propagate through the system leav- 
ing it at contact a with either transmission or reflection at 
any point-like scatterer. Thus C^l represents a classification 
of paths according to contacts and the number of reflections. 
We introduce a partial amplitude (t, E) which describes 
a process during which a particle with energy E enters the 
system through the lead /3, undergoes 2q — 8 a p (for q 7^ 0) 
reflections, and leaves the system through the lead a at a time 
moment t, see Fig.[T] Carriers take a time r = d/v to propa- 
gate from one barrier to the other. This amplitude is a prod- 
uct of amplitudes corresponding to an instantaneous interac- 
tion (reflection/transmission) with time-dependent point-like 
barriers and amplitudes 

i[kd-{e/h) sl j _ T dt'U(t>)\ 

e j (Z) 

corresponding to a free propagation (starting at time tj — r 
and having a duration of r = d/v) between the two bar- 
riers. The time moments, when the instantaneous reflec- 
tion/transmission amplitudes are calculated, are counted 
backward along the path in a descending order starting from 
the time moment t when the particle leaves the system. The 
detailed calculation is presented in Appendix [X] We obtain 



S F (E m ,E n )S F (E n ,E) =I6 m , . (4b) 

n= — 00 

Here 7 is a unit matrix. Substituting Eq. ((1} into Eq. (|4a[) 
we get the following unitarity condition for the matrix 

Sin(t, E)£— 

T 

J^Sl(E,t)S in (t,E)=I. (5) 


Note in Ref.[H3 only the low frequency, r2 — > 0, asymptotics 
for Sin(t,E) was calculated. In contrast equation © is valid 
at any frequency f2 <C /1 of a drive. In the Appendix [B] we 
make a connections between these two calculations. 



III. DC GENERATED CURRENT 

If the Floquet scattering matrix Sf is known then the dc 
current I a flowing out of the structure through the lead a 
can be calculated. To be definite we consider the usual set- 
up for a mesoscopic electron quantum pump. The pump 
connects N r reservoirs. The reservoirs are in stationary 
equilibrium states and have the same chemical potentials, 
fi a — fi, and temperatures, T a — T. Then the distribution 
function f a of reservoir a is the Fermi distribution function 

f a (E) = f (E) = [l + exp (f^)] , where k B is the Boltz- 

mann constant. The current 

f . N r OO 

dEj2 E \SFMEn,E)\ 2 [f p (E)-f a (E n )] . (6) 
•> p=m=—oo 

Using Eqs. CD, © and |A13|l one can calculate the dc current 
generated by the scatterer under consideration. 

To proceed analytically we assume the temperature to be 
smaller than the Fermi energy, 

k B T<^n. (7) 



S in , ag (t,E) e 2i «°f> kd sfr\t,E), (3) 

9 = 

where 2q a/ 3 = 2q + 1 — S a p, k is the wave number, and the 
quantities T (t,E) are given in Eq. (|A13a[) , 

The matrix Si„(t) is periodic in time. This is due to the 
time-periodicity of the frozen scattering matrices L, R which 
describe transmission and reflection at the left and right bar- 
riers. Here we take these barriers to be point-like scatter- 
ers. In addition Si n is periodic in the driving frequency il. 
This follows from the observation that Eqs. ©, (|A13[) give the 
same Si n for r being a multiple of a driving period: r = nT, 
where n is an integer. Therefore, Si n is the same for any 

n = 2-kut- 1 . 

The Floquet scattering matrix is unitary: 

00 

£ S F (E m ,E n )SF(E n ,E)=i5 m ,o, (4a) 



Then calculating Eq. (© we can safely neglect the dependence 
on energy of the scattering matrices of the point-like scatter- 
ers, see Eqs. (|A13|1 . and evaluate them at the Fermi energy, 
E = (i. In addition we approximate the phase factor with a 
linear in energy expansion away from its value at the Fermi 
energy 

e iqkd ^ e iqk F d ^iq—^H-rg ^ 

Thus from now on the propagation time is taken r = To = 
d/vF, with tip the velocity of an electron with the Fermi en- 
ergy. 

With these approximations we perform an energy integra- 
tion in Eq. ((6} and find the dc current as a sum of a contri- 
bution la^ diagonal in dynamical scattering channels and a 
contribution la arising from interference of dynamical scat- 
tering channels, 

I a = /< 0) + ■ (9a) 
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The diagonal part consists of contributions of different scat- 
tering channels that can be labeled by an integer number q 
which is related to the number of times a carrier is reflected 
inside the structure, 

oo 

/i 0) = £4 ?) , Ob) 

<2=0 

with the contribution of the q — th scattering channel given 
by 

= (^(^j^M) . ( 9c) 

The interference contribution is a sum of thermally smeared 
non-diagonal currents 

oo q—1 

/<*> = Re^ J2 e 2l{q - q ' )kFd v(q ~ M'^ , (9d) 

9=1 q'=0 

with 

x 

7(9.9') _ r 



V /act 
/ /\ 2Tvk B T(q-q') ■ . -1 / 27Tfc B Tr /l \ 

V{q-Q) = r - 1 smh B 1 q - q. . 

The factor n describes the effect of averaging over different 
energies of incoming electrons within a relevant energy win- 
dow near the Fermi energy. Note from Eq. (|C2a)l it follows 
that the quantity I a °\ Eq. ()9b|) , is real. 

The time of flight to = <1/vf plays a twofold role in the 
problem under consideration. On one hand, this time sep- 
arates adiabatic T ;§> to and non-adiabatic T < tq regimes. 
On the other hand it defines the energy scale T* = A/(27r 2 fc_e) 
separating low and high temperature regimes. Here we have 
introduced the level spacing A = ttTitq 1 (near the Fermi en- 
ergy) for a double-barrier structure detached from the leads. 
At low temperatures, T < T*, we have 77 = 1, while at high 
temperatures, T>T*, the factor r\ is exponentially small, 

r,*2\ q - q ,\^e-^-«' 1 (10) 

and describes the exponential suppression of the interference 
terms at large temperatures. 

Note the effect of the temperature we consider is only due 
to energy averaging. We do not consider dephasing but as- 
sume that electrons preserve phase coherence as they prop- 
agate through the sample. The crossover temperature T* is 
familiar from the persistent current problem (see, e.g., Ref. l68Tl 
and it appears as a characteristic temperature^ for interfer- 
ence phenomena in single channel ballistic structures. 

The two parts, I a and I a , of the pumped current result 
from different processes that lead to different temperature 



dependencies.— The first part, ft , is a sum of contribution, 
' ar i sm g from different electron's paths inside the system. 
The paths differ by incoming (/3) and outgoing (a) leads, and 
by the index q counting the number of reflections inside the 
system. According to Ref.[67| the quantum pump effect (i.e., 
the effect which results in a dc current generated by the pe- 
riodically driven mesoscopic system) is due to interference of 
various photon-assisted amplitudes describing the interaction 
of propagating electrons with a scatterer that varies period- 
ically in time. Therefore, one can consider the contribution 
as due t° photon-assisted interference processes taking 
place at the same spatial path C^. Each such path can be 
characterized by a delay time 2q a pT, i.e., the difference of 
times when an electron enters and leaves the system. If this 
time is not small compared with the driving period then dy- 
namical effects become important for an electron scattering 
off the system. Therefore, one can consider the path as 
an effective dynamical scattering channel. Then we interpret 
J a as arising due to intra-channel photon-assisted interfer- 
ence processes. Since all the quantum-mechanical amplitudes 
corresponding to such processes are multiplied by same dy- 
namical factor e 2l<3cv ' 3fcd , the corresponding probability is in- 
dependent of energy. Consequently the energy integration 
becomes trivial. 

In contrast, the second part, I a , due to interference be- 
tween different paths (i.e., due to inter-channel interference) is 
defined as a sum of terms oscillating in energy. Consequently 
it vanishes at high temperatures. 

Equation @ gives the dc current generated by the double- 
barrier structure at arbitrary temperature and at arbitrary 
driving frequency. At slow driving, T 3> r, it agrees with the 
famous Brouwer's result 2 for an adiabatic pumped current. 
At larger driving frequencies, it reproduces the numerical re- 
sults of Ref.|7| obtained for a double-barrier structure in the 
zero temperature limit. 

From Eq. ((9| it follows that with increasing temperature 
or driving frequency the different dynamical scattering chan- 
nels contribute independently to the dc current, I a w I a . 
With regard to the frequency this follows from the observa- 
tion that the ratio Ia ^ //« behaves as Otq. Therefore, at 
O — > 00 the contribution dominates. So we can con- 
clude: (i) At high temperatures, T 2> T* the dc current (both 
adiabatic and non-adiabatic) is independent of the Fermi en- 
ergy (over a scale much smaller than the Fermi energy itself), 
(ii) At high driving frequencies, Q ^> tvtq 1 , the dc current 
becomes independent of the temperature. The latter is il- 
lustrated in Fig. [5] (compare with Fig. 3 in Ref.Q), where we 
give the high-temperature dc current I a ' for a double-barrier 
structure with and without the uniform oscillating potential 
U(t) as a function of the driving frequency Q. One can see, 
that the presence of an oscillating uniform potential does not 
change the universal high-frequency behavior of the resulting 
dc current. 



IV. ZERO FREQUENCY NOISE POWER 

We consider the zero- frequency Fourier coefficient of the 
correlation function of currents generated by the dynamic 
scatterer: 
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FIG. 2: (color online) The charge Q = 27r/ (0) /n in units of 
e pumped for a cycle as a function of fi. The black solid 
line is for eU = 0. The red dash-dotted (blue dashed) line 
is for eU = 5 A and ipu = vr/4 (99(7 = 0). The parameters 
are: L = 200tt; V i = V02 = 20; V11 = Via = 10; (pi = 0; 
<^>2 = 7r/2. We use the units: 2m = 7i = e = 1. The Fermi 
energy is fj, = 17.4 . 



Thermal noise 



The thermal noise can be further separated into two con- 



tributions P„' 

ct{3 



,(th,0) 

a/3 



+ p l*,3' !) ' where 



>(tfc,0) _ 2^ 

a/3 

T 



f-^sT (2(5 a/ 3 — (?a/3 — (J/3a) 



g a p = S$gap{t), g afi (t) = £ sj^fo /*) I 



(14a) 



? =0 1 



is the contribution that survives at high temperatures. Note 
that the (dimensionless) matrix g and hence pf" 1 ' ' is inde- 
pendent of the oscillating potential U (t) of the well. 

At low temperatures, T < T* there exists an additional 
contribution P^' 8 ' to thermal noise: 



•y(th,i) 4 



9-1 



nf,2i(q-q')k F d 



i^fcsTReE E [q-q'W 

9=1 q'=0 



1 7=1 



(14b) 



1 OO 

\f> = Jj: J dt'{AL(t)Ai g (t')+Aif}(t')Ai a (t)), (11) 



where A7 = /—(/), and (• ■ • ) denotes quantum-statistical 
averaging. This quantity can be represented as a sum of a 
thermal (Nyquist- Johnson) noise P„a' an d a shot noise P„g , 

Pa/3 = P^a + Pa/3 • Notice, each of these contributions 
depends on both the temperature and the drive. The reason 
to divide the whole noise into two parts is rather conventional: 
The former contribution vanishes at zero temperature while 
the latter one vanished in the absence of a drive. In terms of 
the elements the matrix Sin it, E) the noise is^i 



_ 2e z k B T 

a/3 ~ 



V 



I din, 0a 

(t, E)\ 2 + 6 a p |l + Ein,aa{t, t] 



(12a) 



^ = e iSdE £ [fo(E)-ME m )fff^ 



m— — 00 



(12b) 



xe"" 1 ' 2 ll) S* nMf3 (ti,t2\ E)S irlia p(ti,t2; E m ) ■ 
Here Sin is the two-particle scattering matrix^ 

E m {ti,t 2 \E) = S m {ti,E)S\ n {E,t 2 )- 

To calculate the noise power generated by a driven double- 
barrier structure we use Eqs. J3J, ©, (H} and (IA13|) . 



(13) 



The function r)(q — q') is defined in Eq. (|9cl) . 

At high temperatures the thermal noise of a driven sys- 
tem depends on the matrix g in the same manner as the 
thermal noise of a stationary scatterer depends on the con- 
ductance matrix— Therefore, we can treat g as an effective 
(high-temperature) conductance matrix of a driven scatterer. 

At adiabatic driving, T 2> to, the matrix g^ ad ^(t) = 
Km g(t) is independent of Q. It can be expressed in terms of 

reflection/transmission coefficients for point-like scatterers. If 
L and R are symmetric in lead indices then we get, 



1 



gif (t) T L {t) Tnit) 



1 . 



(15) 



Here T L 



Tr = \R 12 \ 



and g^a \t) 



Note that the effective transmission coefficient g[ 2 d \t) dif- 
fers from the probability for incoherent (sequential) tunneling 

through two barriers, g[ 3 2 q) it) = [T£ l it) + T' 1 ^)] ~\ The 
latter is attributed to structures with inelastic relaxation pro- 
cesses inside the well^ whereas Eq. (|15[) can be obtained as a 
consequence of incoherent but energy conserving scattering 74 
in the well. Note also, that in Ref.[7a it was shown that in the 
adiabatic regime the dc current pumped by the double barrier 
structure with strong inelastic scattering is determined by the 
sequential conductance (Ji^ 9 '- l n contrast the high tempera- 
ture pumped current Eq. ()9b|l . at SI — > is determined 
by the quasi-elastic adiabatic conductance ff^ • 

At high driving frequencies g becomes dependent on Q. In 
Fig. [3] we give the transmission probability g± 2 as a function 
of the non-adiabaticity parameter Qtq. 
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At high temperatures, T » T" 




FIG. 3: The effective transmission probability g±2 as a func- 
tion of fl. The parameters are the same as in Fig. [2] 



B. Shot noise 

Similar to the thermal noise, the shot noise has a part which 
vanishes at T ^> T*. To show this we use Eq. <(3j and repre- 
sent the two-particle scattering matrix Ei n , Eq. (|13[1 . as fol- 
lows: 



Sn,a^(ti,* 2 ;B) = E E e 2iQ «V d r« (ti,t 2 ;£), 



(16) 



Substituting this equation into Eq. (|12b|) and integrating 
over energy we obtain the following: 



Pi^ = ¥ E E e 2ik * dQ $^V% \ qi ,q 2 ) , 

91 ,92 = -oo 7,<5 

<p(sh) _ jj dtidt 2 e *"i(i2-tl+fim©^ > [ 91 '' 32 (I'D 



m=-oo 



where 

F% 92) = V (e^ 3 [«i,«.]) { - 2fc s Tcos (mftn,e^t9l.«>: 



+ft coth 



\ m ^ B T) T0 e^ lqi , q2] / 



contribute to the shot noise. In this limit, for 



instance, the auto-correlations are 



only terms with 



,(sfi,0) _ 2e 2 



rft i dt2 

t 

9=0 



TrfJ^ E iMi 2 -ti + 2gfiro) 



x(2-«5,o) 



E ^aan (^1 1 ^ 2 > A 4 ) 
7 



F T (i) = 2fc B T E cos(mi) coth (g™) - l} 

771—1 \ / J 



(18) 



From Eq. (JTSj it follows that for ftfi < fc s T the shot 
noise is inversely proportional to the temperature: p( 3 ' 1 ' ' ~ 
(fcsT 1 ) -1 . The dependence on frequency il is more subtle. 
For adiabatic driving 17 <C Tq" 1 only Ft depends on the 
driving frequency, Ft ~ fi 2 . Therefore, in this case we get 
p(sfi,o) ^ Q?/(k B T) in accordance with Refs. [rOEl. In con- 
trast, for non-adiabatic driving, t^ 1 <C <C ksT/h, the 
matrix Sin and hence E oscillates in f2. As a consequence the 
shot noise oscillates in fi with increasing (~ fi 2 ) amplitude. 
Since HQ <C ksT the thermal noise dominates over the shot 
noise in both the adiabatic and the non-adiabatic regimes. 

The situation is different in the strongly non-adiabatic 
regime, r -1 <C ksT <C hfl. In this case Ft ~ fi. Therefore, 
the shot noise is temperature independent and it oscillates 
in frequency with amplitude growing linearly in SI. In this 
regime the shot noise exceeds the thermal noise. 



V. THE EFFECT OF AN EXTERNAL BIAS 

In this section we extend the formalism presented above to 
the case when there is in addition an external bias. Several 
factors make such a consideration important i£S First, usually 
the mesoscopic structure is a part of some electrical circuit 
with non-zero impedance used to measure a current. In the 
presence of a pump current the voltage arising at the exter- 
nal impedance acts back on the sample. Second, the gates 
with oscillating voltages used to drive the sample can produce 
oscillating voltages at all the terminals.— The rectification 
of ac currents due to external ac voltages by the oscillating 
scattere r 19 ' 77,78 ' 79 can lead to a dc current in addition to a 
pumped current. 



A. The effect of a dc bias 

To be definite we assume that the dc voltage eV is applied 
to the right reservoir: 

Hi=H, fi2 = H + eV, eV<^fi. (19) 
We calculate the dc current, 7i, via Eq. ©. In the high- 
temperature limit, I » T* we find: h = lf ] + l[ dc) ■ Here 
Jj ' is the pumped current, Eq. (|9a[) , and I^' is a current 
due to a dc voltage: 



l[ dc) = VG Q gi2 . 



T > T* . 



(20) 



where Go = e 2 /h is the conductance quantum (per single spin 
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channel), the effective conductance g\2 is defined in Eq. (|14a|l . 
Note that the pumped current and the dc volt age- driven cur- 



rent are conserved separately, E, 



r(0) 
1 a 



and J2 a I, 



(dc) 



0. 



In the two terminal case of interest here, the last equation 
leads to gi2 = f/2i- In addition using Eqs. (|14a|l and (|C2a|l we 
find that the sum over all the incoming leads equals to unity, 
E/3 9<*p{t) — 1- This condition allows us to interpret g a p{t) 
as equal-arrival-time transmission (for a 7^ 0) or reflection 
(for a = jS) probabilities for the particles with Fermi energy. 
In other words, the quantity g a 0(t) in Eq. (|14a|l . is a probabil- 
ity for an electron entering the system through the lead f3 to 
leave the system at a given time moment t through the lead a. 

I for 1 I 2 

This probability is a sum of probabilities, p,Tg (t, fi)\ , cor- 
responding to electrons entering the system at time moments 
t — To(2g + 1 — 8 a p)- We emphasize that this interpretation 
applies at high temperatures, T 2> T* . At lower tempera- 
tures only the relation Eq. ((5j), averaged over a pumping cycle, 
holds. 

In the non-adiabatic, high-temperature regime, Qto 3> 1 
and T 3> T* , both the pumped current and the dc 
voltage-driven current oscillate at the driving frequency 
(see Fig. [2] and Fig. [3} with the same period. The important 
difference is that the amplitude of oscillations of grows lin- 
early with Q. The pumped current can be expressed in terms 
of the average charge transferred per cycle Q(Q.) multiplied 
by the frequency, i' ' = QCl/(2ir). The charge Q pumped 
per cycle as a function of the pump frequency Q is given in 
Fig-El The amplitude of oscillations of /' dc ' is linear in the 
dc voltage V. Therefore, varying both f2 and V one can, in 
principle, distinguish between the pumped current and the 
external voltage-driven current. 

We reemphasize that the quantities g a p define both the 
thermal noise and the dc conductance of a driven scatterer. 
In contrast the rectification of currents due to external ac 
voltages is defined by different quantities. 



B. The effect of an ac bias 

Let us suppose that ac voltages are applied to the contacts 
of a dynamic scatterer (a — 1,2, ■ ■ ■ , N r ): 

V a {t) = Vo, a + V a cos(Q,t + <pv a ), eVo, a , eV a -C fMc, ■ (21) 
Then the dc current flowing through lead a is^ 



I a = %fdEi -/«(£)+£ E ME-nnn) 

( /9 n=-oo 



x E Sp ia 0(E,E q )SF }a f3(E,E rn )Yp n+q Y l 3 >n+m 

m,q— — 00 



(22) 



Here Tp,n is the Fourier coefficient for the function Tp(t). 
The latter is given by Eq. ()A4|) with ell and ipu being re- 
placed by eVa and (pvg , respectively. We use the symbol tilde 
to indicate that the dc current I a is calculated in the pres- 
ence of external ac voltages. The Fermi functions f a (3) enter- 
ing Eq. (|22[) depend on the corresponding chemical potentials 

Ma(/3) = M + eVo, a (/3)- 



To simplify the equation given above, we introduce the ma- 
trix Sout(E,t) whose Fourier coefficients define the elements 
of the Floquet scattering matrix in the following way>££ 



S ou t,n(E) = Sf(E, E- n ) . 



(23) 



The lower index "out" indicates that the matrix introduced 
in this way depends on an outgoing energy entering the cor- 
responding element of the Floquet scattering matrix. In con- 
trast the matrix Sin which we introduced earlier is a function 
of an incoming energy. 

Substituting Eq. (|23l) into Eq. (|4b[) we derive the following 
constraint to the matrix S ut{E,t) [compare with Eq. Q]: 



1 

! 



-^:Sout(E,t)Sl ut (t,E) = 7. 



(24) 



To calculate this scattering matrix we find Sf{E„, E) using 
Eqs. (TTJ, © and l(AT3]l . Then with the help of Eqs. (T53J after 
the inverse Fourier transformation we find: 



S out . a p{E,t) = £ e 2 ^ kd S% T \E,t) . 

q = 



si q f\E,t)=sly ) (t + 2g cl pT,E). 



(25) 



The amplitudes S^f\t,E) are defined in Eq. (|AT3|l . We in- 
terpret S ou t >a p{E,t) as an amplitude corresponding to scat- 
tering of an electron which enters the system through the lead 
j3 at time moment t and leaves the system through the lead a 
with energy E. Since the energy of the outgoing state is fixed, 
the time when the particle leaves the system is not fixed. Cor- 
respondingly, since the time moment when the particle enters 
the system is fixed, the incoming energy is not fixed. 

Now we proceed with the calculation of the current. After 
simple transformations we rewrite Eq. (I22II in the following 
way: 



r —1 + r' 



(ex) 



(26) 



where the current I a is the pumped current, Eq. in the 
absence of the oscillating voltages at contacts. The excess 



current Ic\ x ^ due to ac voltages is 



/r^t/d^E e [/„(£_„)-/«(£)] 

j3 n— — 00 

x jj(*^out,a/3^/3) n | — \S ut,af3,7i\ j" • 



(27) 



The effect of oscillating potentials at external reservoirs is 
formally only a phase shift [see Eq. (|A4|| ] of the correspond- 
ing elements of the scattering matrix Sout- However due 
to the directional asymmetry of time-dependent scattering, 
Sou*, a/3 7^ S out , /3a, the current I^^ can be non-zero even 
if the voltages at all the contacts are the same, V a (t) = 
V(t). Therefore, in accordance with the conclusion reached in 
Ref.HH the excess current contains two generic contributions, 



r(rect) - r(int) 
l a -y- 1 a 



The first one, /, 



(rect) 



results from 



rectification of ac currents due to external ac voltages by the 



oscillating scatterer. The second one, , is due to inter- 

ference between internal and external ac currents. On the 
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FIG. 4: (color online) The transmission probabilities c/2i(t) 
(black solid line) and §12 (i) (blue dashed line) are given 
as a function of time for one pumping period. The non- 
adiabaticity parameter is Qto = 7T. Other parameters are 
the same as in Fig. [2] 



FIG. 5: (color online) The excess I{ (black solid line), rec- 
tified /( rect ' (blue dashed line), and interference, J^ nt ^ ( re d 
dash-dotted line) currents in units of Jq = VGo as a func- 
tion of fi. The external voltages are V2 = V cos(Q,t), V\ = 0. 
Other parameters are the same as in Fig. [2] 



other hand one can think of the latter contribution as due to 
external voltages which act as additional pump parameters. 
Note that the interference contribution was also addressed in 
Ref-ETj. 

In the high-temperature limit Eq. (|27|l can be simplified. 
Using Eq. (|25|l and performing integration over energy in the 
same way as we did to get Eq. © , we obtain: 

T 

lL ex) = Go f ^ Yl Mi)9<* {t), T » T* , (28a) 



~g aB (t) = Y,\s% T0 \t^)\ ■ (28b) 

9 = 

The pumped current I a , Eq. ©, is conserved: ~}2 a I a = 0. 
Therefore, the conservation of the whole current Eq. (|26|l . 

implies 1^ = 0. The latter in turn leads to the fol- 
lowing constraint: g a p(t) — const. Quite similar to the 
case with g a fj we find that this constant is the unity. Using 
this condition one can split 7„ ' into the rectification [i.e., 
vanishing at Vi(t) = Va(t)] and interference parts as follows: 

= Go J f [v 2 (t) Vl (t)] ^+92^) (29a) 


r 

r (int) n f dt Vl(t) + Vl(t) r_ , . . .-, . . 

1} ' = Go / -- 2 b 12 W ~ 521 wJ • ( 29b ) 



The interference contribution can be non-zero, j£ m *' ^ 0, 
if and only if the dynamical scatterer shows different time- 
dependent transmissions to the left and to the right, g\2 (t) 7^ 
521 (t), see Fig. [4] In addition, the time-dependent voltages 
should be present, V a (t) ^ const. The last statement follows 



from the fact that after averaging over time the matrix g 
becomes symmetric in lead indices. Therefore, we conclude 
that the current /' lnt ' j s due to an interplay of the external 
induced dynamics and the internally induced dynamics. 

The excess current I a ex \ arising in the presence of exter- 
nal ac voltages, comprises two contributions: a pure recti- 
fication current, I a rect \ and an additional pumped current, 
Ia ■ These two contributions differ from each other. In ad- 
dition they differ from the current pumped by the unbiased 
dynamical scatterer. Formally, the last mentioned difference 
is due to the fact that is defined by the scattering ma- 
trix Sout while I a is defined by the matrix Sin- This results 
in a different frequency dependence, see Fig.[5] for and 
Fig. [2] (black solid line) for the pumped current. Everything 
mentioned makes it difficult to extract a pure pumped current 
from an experimentally measured dc current. 

If the voltages at contacts are time- independent, Vp(t) = 
const, equation (|28a|l transforms into Eq. (|20|l since the av- 
erage over a pump period for quantities g a p and g a p is the 
same. This follows directly from Eq. (|25[) (the second line). 

Note by analogy with g a p(t) we interpret g a p(t) as equal- 
departure-time probabilities for transmission/reflection of 
electrons with Fermi energy: This quantity is the probabil- 
ity for an electron entering the system at a given [defined by 
the potential Vp(t) in Eq. (|28a[l ] time moment t through the 
lead 13 to leave the system through the lead a. Accordingly 
to Eq. (|28b|) the events with electrons leaving the system at 
different time moments do contribute to g a g(t)- We note the 
equation (|28[) was obtained for high temperatures. For lower 
temperature we have no instant-time quantity like g a /s (t) with 
a transparent physical meaning. The matrix Sout entering 
Eq. (|27|l obeys a constraint that is non-local in time, Eq. (|24|) . 
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VI. DISCUSSION AND CONCLUSION 

We have analyzed current and noise generated by a dy- 
namical double-barrier structure with and without an exter- 
nal bias at arbitrary temperature and at arbitrary frequency 
and amplitude of a periodic harmonic drive. 

To perform such an analysis we have developed a method 
to calculate analytically the scattering matrices Si n (t,E) and 
S ut(E,t) whose Fourier coefficients define the elements of 
the Floquet scattering matrix, Eqs. fT) and (|23[l . For the 
elements of these matrices which we found for the double- 
barrier structure there exists a simple interpretation. This 
interpretation allows to find the corresponding matrices for 
more involved systems. For instance, the matrix element 
S in , a/3 (t,E} = Y. q e 2l " c ' 0kd S% T \t,E), Eq.©, can be calcu- 
lated as a sum over all the possible paths which an electron 
can follow to propagating through the system. Here the par- 
tial amplitude e 2iqa i 3kd S^ l f^ (t, E) corresponds to a particle 
entering the system with energy E through the lead /? and 
leaving the system through the lead a at a time moment t. 
The summation index q corresponds to the path C ^l having 
definite length (in our case for q > this length is a product 
of the number 2q — 5 a/ s of reflections inside the system and 
the distance d between the two barriers). Each partial ampli- 
tude is a product of amplitudes corresponding to an instanta- 
neous interaction with time-dependent point-like barriers (in 
our case these amplitudes are the elements of matrices L and 

R) and a number of amplitudes e^ kd eh t i~' r dt U ^ " corre- 
sponding to free propagation (having a duration of r = d/v) 
between the barriers. The time moments, at which the cor- 
responding amplitudes are calculated, are counted backwards 
along the path in a descending order starting from the time 
moment t when the particle leaves the system. 

We found an analogous interpretation for elements of the 
matrix S ou t(E,t). The in- and out- scattering matrices be- 
come equal (and equal to the frozen scattering matrix) in the 
limit of a slow drive, f2 — > 0. At any finite driving frequencies 
(£2 ^= 0) they are different. Nevertheless there is a simple re- 
lation between these matrices. To find it we use the micro re- 
versibility of the equations of motion. Let the driving parame- 
ters depend on time in the following way: Pi(t) ~ cos(fit+¥?i), 
see, Eq. ()A1 [) . Then the elements of the Floquet scattering 
matrix are subject to the following symmetry^ 

SfME, E n ; H, ip) = S Fl0a (E n , E- -H, - v ) . (30) 

Here tp is the set of all the ipi's. In addition we introduced a 
possibly present magnetic field H . Substituting Eqs. (|A9[) and 
(|23|) into Eq. (130[) and performing the inverse Fourier trans- 
formation we arrive at the desired relation: 

Sin t0 !p(t, E; H, ip) — S ou t t p a (E, —t; —H, —tp) . (31) 

Note that each term of the decompositions, Eq. © and 
Eq. (f25)l . satisfies Eq. (pi) separately. 

We showed that the scattering matrix Si n defines currents 
generated by the oscillating scatterer itself, while the dual 
scattering matrix Sout defines currents flowing through the 
scatterer under an external bias. Here we calculated the 
matrix Sin for arbitrary frequency and consequently we can 



calculate the current in the adiabatic as well as in the non- 
adiabatic regime. 

The adiabatic regime can be defined^ by considering the en- 
ergy dependence of the stationary scattering matrix . Since 
both transmission and reflection times are of the order of the 
dwell time, the latter is (at least in the absence of interac- 
tions) the relevant time scale. The results obtained here al- 
lows us to conclude that the dwell time is appropriate to de- 
fine the adiabatic regime for zero-temperature pumping. In 
contrast, at high temperatures another time, the propagation 
time to = d/v (necessary for an electron to travel between the 
two barriers separated by the distance d), defines a relevant 
time scale. 

Let us consider the scattering processes which are relevant 
for pumping. Generally speaking the interplay between the 
quantized energy exchange and the interference of photon- 
assisted scattering amplitudes is at the origin of the pump 
effect.^ As we showed, there are two types of interference 
processes. The first one takes place within the same spatial 
path (intra-channel interference). A second process includes 
the interference between different spatial paths (inter-channel 
interference) . The former contributes to the generated current 
at low as well as at high temperatures, while the latter con- 
tributes only to the low temperature pump effect. The dwell 
time (defined in the stationary case) depends essentially on 
interference processes of the second type (there it nothing to 
interfere within the same spatial path in the stationary case) . 
Therefore, we expect that the dwell time is relevant only to 
low temperature pumping. In contrast for high temperatures 
we expect that a time of the order of the propagation time to 
defines the lower bound for the period of a drive separating 
adiabatic and non-adiabatic dynamical regimes. 

Our example suggest that distinction of inter- and intra- 
channel interference processes is also important for the mag- 
netic field symmetry of the current pumped by a two-terminal 
scatterer. The intra-channel interference processes are not af- 
fected by the magnetic field. However the inter-channel pro- 
cesses are sensitive to the magnetic field. Therefore, we can 
conjecture that with increasing frequency and/or tempera- 
ture (when only intra-channel interference processes matter) 
the pumped current has to become symmetric in magnetic 
field like the stationary (two-terminal) conductance. This is 
in striking contrast with the low temperature adiabatic pump 
effect which has no definite symmetry under the magnetic field 
reversal, see, e.g., Refs. 00- Our conjecture agrees with ex- 
perimental results presented in Refs. [l9ll2ll for large pumping 
frequencies and powers. 

The method presented in the present paper is also useful 
to calculate the time-dependent current generated by the dy- 
namical scatterer. Actually it was used in Ref.H^ to describe 
the nonlinear response of a coherent capacitor to a large am- 
plitude drive. In a recent experiment^ it was demonstrated 
that such a capacitor can serve as a very fast single electron 
source. 

In conclusion, we developed a method to calculate ana- 
lytically the Floquet scattering matrix of a dynamical meso- 
scopic system consisting of a set of point-like scatterers and 
uniform potentials in between them. Our approach allows to 
go beyond the adiabatic approximation and to analyze the 
properties of the generated current and its noise at arbitrary 
frequency and amplitude of the drive. 
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To match the wave function coefficients at different spatial 
regions we use the boundary conditions formulated in terms 
of the scattering matrices Lf and Rf for point-like scatterers 
located at x — and x — d, respectively. These conditions 



APPENDIX A: THE FLOQUET SCATTERING 
MATRIX 

Let us consider a one-dimensional structure consisting of 
two oscillating 5-function barriers and an oscillating spatially 
uniform potential in between the barriers. All potentials os- 
cillate with the same frequency. The single-particle scattering 
at a periodically driven system is described via the Floquet 
scattering matrix Sf^— To find its elements one needs to 
solve the Schrodinger equation for an electron wave function 
^E{t,x) with, in our case, the following time-dependent po- 
tential: 



V(t,x) = V L (t)S(x) + V R {t)8(d-x) 

+eU(t)0(x)6(d - x) , 
V L /n{t) = V 0L/R + 2V 1L/R cos(Q,t + tpi/ R ) , 
U(t) = U cos(«t + ipu) , 



(Al) 



where 9{x) is the Heaviside step function. 

To calculate Sf,u and Sf,2i we consider a plain wave of a 
unit amplitude with energy E incident from the left. Then the 
wave function ^e{P,x) outside the scatterer can be written 
as follows: 



y E (t,x) = e-'n t J2 e 

n — — oo 



1pn{x) = < 



' 5n ° etkX + \[i S pM E n,E)e- ihnX ,x<0, ( A2 ) 



Here E n = E + nhSl, k n = \/2m e E n /h, with m e being the 
electron mass. 

Following Ref.|65| we represent the functions tp n {x) inside 
the scatterer in the form: 



oo 

M*)= T «" ! (ne klX + bie~' k ' x ) , 0<x<d, (A3) 



where the coefficients ai, bi are independent of x and t. Here 
T n is the Fourier coefficient for a function T(t) describing the 
effect of an oscillating uniform potential U(t), 



T n = J f e innt ?(t) , 



T(t) = exp [-ifg sin(m + <p v )] , 
where T = 2n/Q is the period of the drive. 



(A4) 



>£SFMEn,E) = ^L F ,ll(En,E) 

oo j oo 

+ X] J^-iF,12(-En, E m ) T m _;6; . 

m— — oo / — — oo 



^2 T n - P a p = <J -^Lp,2i{E n ,E) 



+ E J^L F MEn,E m ) E r m -ibi. 

m— — oo * I — — oo 



(A5) 



p— — oo 



Yl \ jf^^-F^niEn, Em) T m -iaie 



i k j d 



m— — oo 



T-SF,2l(En, E) 



E 



m— — oo 



We suppose that the frequency fi of the drive is small com- 
pared with the relevant electron energy E, 



no. < E. 



(A6) 



We solve Eq. (|A5|I to zeroth order in hQ/E. Within this accu- 
racy we can simplify Eq. (|A5[) . First, the Floquet scattering 
matrix Xf(E„, E) for a point-like scatterer (i.e., for a scat- 
terer with a spatial extend much smaller than the de Broglie 
wave length for an electron with energy E) can be expressed in 
terms of the Fourier coefficients X n (E) for a frozen scattering 
matrix X(t, E) = X(t + T, E) (i.e., the stationary scattering 
matrix with a strength being dependent on time)>£2 



Xf{E u , E m ) = X n - m (E) . 



X = L,R. 



(A7) 



For a single 5-function barrier V(t,x) = V(t)S(x) the frozen 
scattering matrix is well known. Its elements are X a/ s(t, E) = 
k/[k + iV{t)m e /h 2 ]-5 a p. 
Second, we can put 



fan 
k' ri 



ikd HVlT 

e e 



(A8) 



where r = d/v is the time necessary for an electron with 
energy E to propagate from one barrier to the other barrier. 

Next we use that the quantities a;'s and fei's are Fourier 
coefficients for some functions a(t) = a(t + T) and b(t) = 
b(t + T), periodic in time. In addition we introduce a matrix 
Si„(t,E) periodic in time whose Fourier coefficients Si n ,n{E) 
define the Floquet scattering matrix Sp(E n , E) for the whole 
structure as follows: 

S in ,n(E) = S F (E n ,E). (A9) 
The lower index "in" indicates that this matrix is a function 
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of an incoming energy. With these definitions one can easily 
perform the inverse Fourier transformation of Eqs. (|A5|l and 
arrive at 



S in ,i 1 (t,E)=L 11 (t,E) + L 12 (t,E)T(t)b(t), (AlOa) 



T{t)a(t) = L 21 (t,E) + L22(t,E)T(t)b(t), (AlOb) 



„(8>0) 



q-1 



(*) = II Ru{t 2] )L 22 {t 2j - r) 



3=0 



(A14) 



= 1 . 



In the next section we consider in detail adiabatic (to first 
order in Q) asymptotics for the scattering matrix derived 
above. 



' l T{t)b{t + t) = e ikd R n (t, E)T(t)a(t - r) , (AlOc) 



S in ,2i (t, E) = e lkd R 2 i (t, E)T(t)a(t - r) 



(AlOd) 



Substituting Eq. (|A10c|l into Eq. (|A10b|l we find the non- 
local in time equation for the function a(t) (for brevity we 
suppress the argument E): 

a(t) = T*(t)Lai(t) + e 2kd L 2 2(t)Rn(t - r)a(t - 2r) , (All) 



where the star denotes complex conjugation. Since the coef- 
ficients entering the above equation with amplitudes smaller 
than unity, we can formally write down the solution in the 
form of an infinite series, 



a(t) = E e l2qkd X M (t)T*(t - 2qr)L 21 (t - 2qr) . 

9=0 

AC«*»(t)= ULi^Rnitv-r). 

3=0 



(A12) 



Here t 2j = t - 2jr. At q = we put A (0) (t) = 1. 

Substituting Eq. (IA12[) into Eqs. (IA10|l we obtain the ma- 
trix elements Si n ,ai, a = 1,2. To find the matrix elements 
Si„, a2 a problem with a plane wave of a unit amplitude com- 
ing from the right has to be considered. The solution to these 
two problems is written in Eq. ((3j with 



(t,E) = e-^aiY\t,E) 



(A13a) 



where *«f> = eft" 1 f t _ T2q ^ dt'U(t') and 

ffg' T) (t)=Lii(t), 
a[l >0 ' T \t) = L 12 (t)R 11 (t-T)L 21 (t 2q )\^- 1 \t-2r), 

aif T \t) = R 2 i(t)L 21 (t 2q+1 )\ (q Ht~r), (A13c) 



(A13b) 



(t) = L 12 (t)R 12 (t 2q+1 )p (q) (t-r), 



42 T) (t)=R22(t), 

a ( 2 f°- T) (i) = R21 (t)L 22 (t - T)R 12 {t 2q )p^ (t - 2r). 
Here the function p' 9 ' (t) is 



(A13d) 



(A13e) 



APPENDIX B: ADIABATIC DRIVE: fir < 1 

In this Appendix we verify that for a slowly driven scatterer 
the matrix Si n = +0(ft 2 ) can be represented through the 
frozen scattering matrix S as it was proposed in Refs. [52II66I : 



S%(t,E) = S(t,E) + 



ih d 2 S 
2 dtdE 



+ hflA(t, E) . 



(Bl) 



We calculate the matrix A which is responsible for generat- 
ing a directional asymmetry in scattering on an adiabatically 
driven scatterer. In turn, this asymmetry leads to an adia- 
batic quantum pump effect. The matrix A satisfies^, 



Ml ( S^A + A* S 



'P{S' t ;5'} = ih 



dff 1 " dS _ dS f d§ 
dt dE dE dt 



(B2a) 



(B2b) 



This equation results from the current conservation law up to 
linear in Q terms. 

To find the matrix A we use the anzats Eq. (|B1[) and con- 
sider the solution Eq. Q in the limit of low driving frequency 
Q — > 0. In this limit the delay time r is small compared with 
a driving period T = 2tt/uj. We use the ratio e = r/T as a 
small parameter. 



Zeroth order 



In zeroth order in e the matrix coincides with the 
frozen scattering matrix S. The frozen scattering matrix is 
the stationary matrix with parameters being dependent on 
time. To calculate the frozen scattering matrix one can use 
Eqs. ((3J (|A13|l and neglect changes in all the quantities during 
the delay time r. Therefore, we have to put t 2q « t in all the 
terms in Eqs. (|A13[) and use: 

~ eU(t)Th- 1 (2q + 1 - 5 a0 ). 
As a result we find (t, E) = S(t, E) with matrix elements: 



S aP {t,E) = J2 S^(t,E) : 



9=0 



(B3) 



S%{t,E) = e *(*«»-trWr/*)C2«+l-«a/j) ff tal»( ti ^ j 

where the matrix <t (<7 ' t) (t, E) is defined in Eq. (|A13|I . 
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2. First order 

To calculate the matrix Sj^ we expand Eq. (|A13[) in pow- 
ers of r and keep only the terms linear in r. Our aim is to 
calculate the irreducible part of , i.e. the matrix A. To 
this end we use Eq. (|B1|) with the matrix S given by Eq. (|B3|I . 
We take into account that the frozen scattering matrix S de- 
pends on time t through U(t) and the matrices L and R of the 
point-like scatterers. The energy dependence of this matrix 
is due to the phase factor e 2 * 9 only. We neglect the energy 
dependence of the matrices L and R and take them at E — fj,. 
We can do so since L and R change on the energy scale of 
the order of fi while in our problem the much smaller scales 
of order HQ and ksT are relevant, see, Eqs. (|A6f) and 0. 

After a straightforward calculation we find 

oo 

hQA a0 (t, E) = J2 S ( a q) B (t, E)A% (t, ti) , (B4a) 

9=0 

where 



per we use the resulting conditions to simplify the expressions 
for the pumped current. Substituting Eq. @ into Eq. (|4bp 
and performing the inverse Fourier transformation we deduce 
(as before we treat 5"' T ' as energy independent over the scale 
of order HQ) , 



E S (q - T) (t,E)S^' T) \t,E) 

9=0 

OO oo 

+ E E e- 2lakd S {p ' T) {t,E)S {p+3 ' T) \t + 2r8,E) (CI) 

p=0 s=l 

+ E E e 2 " kd S {q+s ' T) {t 1 E)S {q < T) \t~2sT,E) =1. 

9=0 3 = 1 



These identities have to hold for any wave number k. Taking 
into account that in our case the matrix S^ q ' T ' can be kept en- 
ergy independent over the interval corresponding to a change 
of kd by 2-7T, we derive next three equations from Eq. (|C 1 [) : 



A (q) 



A (q) 



ro(2q + l) 



Of 



r (2g + l) d 
2 dt 



In 



In 



r 2 J 2 dt n \ L22 



R12 

L12 



T g d 



In 



L22 
Rn 



.(g) d ( R2I 



(B4c) 



(B4d) 



(B4e 



At U(t) = 0, the above expressions coincide with those ob- 
tained in Ref.l52l 

Equation (|B4|I shows that the matrix A possesses symmetry 
properties with respect to interachange of lead indices which 
are different from those of the frozen scattering matrix^ The 
symmetry of the A matrix depends on differences of the ma- 
trix elements of the L and R matrices [see Eqs. 



(|B4c|l and 

l|B4dfl ]. The main point is that the symmetry of the A- matrix 
is fundamentally different from the (frozen) scattering matrix 
symmetry. 



APPENDIX C: ALTERNATIVE UNITARITY 
CONDITIONS 

We can obtain another (alternative) formulation of the uni- 
tarity condition for the matrix Si„ . In the main part of the pa- 



J2S {q ' T) (t,E)S (q ' T)t (t,E) = i , 



9=0 



J2S (p '^(t,E)S 



(j>+s,T)t 



(t + 2rs,E) =0, 



p=0 



9=0 



(9+s,t), 



{t,E)S (q > T) \t-2rs,E) =6, 



(C2a) 



(C2b) 



(C2c) 



where is a null matrix. Notice that while the unitarity 
condition Eq. ([5]) is quite general, the untitartity conditions 
obtained above are valid only within the approximations we 
made to get the matrix Si n , Eq. ©. 

Taking a time derivative of Eq. (|C2a|l we immediately 
prove that Eq. (|9bp is real. In addition using Eqs. (fC2b) and 
(|C2c|) we can show that the product of scattering amplitudes 
5 (9 ' To) (t,M)5 (9 '' To)t (t - 2r [q - q'],^) (corresponding to par- 
ticles leaving the scatterer at different time moments, t and 
t — 2tq [q — q']) in fact does not contribute into an interference 
W 



current l£> , Eq. (gd 
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